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ABSTRACT 

A baroclinic model programmed for numerical solution is used to in- 
vestigate the effects of diabatic heating on ene development of pressure 
waves for periods up to 24 hours. Several forms of the model are studied 
with results being obtained for various wavelengths. The results lead to 
some meteorological improbabilities; however, the relative development of 
the pressure wave appears valid for short time periods. 

The writer wishes to express her appreciation to Professor George 
Haltiner of the U. S. Naval Postgraduate School for his assistance and 
guidance in this investigation. Appreciation is also expressed to LT 
Jack Nourie, USN, for his assistance in preparing diagrams and compiling 


data. 
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LIST OF SYMBOLS 
Variables and Constants - Roman Alphabet 
amplitudes of various sine and cosine functions 
Rossby wave speed,= ({ — Pf > 
shorthand equivalent of various equations as 
designated in text, (appear with or without primes) 
the Coriolis parameter, QS? sin 
an average value of the Coriolis parameter 
acceleration of gravity 
unit vector in X direction 
unit vector in Z direction 
wave number, 2/7] /L 
wave length | 
atmospheric pressure 
heating term 
amplitude of the pressure wave 
time 
mean zonal flow in X direction 
component of perturbation velocity in X direction 
horizontal velocity vector 
component of velocity in Y direction 
coordinate toward east 
coordinate toward north 


height of isobaric surface 


iv 
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Variables and Constants - Greek Alphabet 


aw 
kA p 
coefficient of certain terms, Date 46 ( £ ) 
Oo 


variation of the Coriolis parameter, Of/d y 
j : = 
vertical component of relative vorticity, mY. Dp 
ao 
phase angle 


potential temperature 


2d 20 


tability factor, =—-+ 
stability factor © 9P 3? 


geopotential, gz 

geographical latitude 

total derivative of pressure with respect to time 
the angular speed of the earth's rotation 


positive tolerance 


Subscripts and Superscripts 
vertical index denoting a pressure level 


index denoting the nth estimate of a quantity 


Operators 
two dimensional del operator on a constant pressure 
surface 
Laplacian operator 
finite difference operator with respect to p and t 


absolute value of a quantity 





l. Introduction 

Most numerical prediction methods for obtaining prognostic charts of 
the pressure field neglect diabatic heating effects under the assumption 
that their influence is negligible for short period prognosis. 

Under certain circumstances, however, the exchange of sensible heat 
between an air mass and the underlying surface may be considerable. It is 
therefore the purpose of this paper to investigate the effects of such 
diabatic exchange on the development of pressure waves for periods up to 


24 hours using a numerical model. 





2. The Model 


The model with which the aforementioned heating effects were studied 
was developed from the geostrophic relationship, the vorticity equation, 


and the thermodynamic equation, 


cE XS, (1) 
oF + V°TV(F +5) = 6 aw , (2) 
Se oP 
ao +V'Vad + 6w=-Q. (3) 
opot oT ei 


The dependent variables are expressed as the sum of a base state plus 


perturbation quantity (indicated by a prime), 


b= Ph Cp,4) +P (%4,7,4), (4) 
\ =e Ce) @ + V (¥,4,7,% ) (5) 
Cee eG Kh, Pe) (6) 
Q= OCP) + Ole) with @¢p) =O, (7) 
$= §+6, (8) 


Equations (2) and (3) are then linearized by the usual method of assum- 
ing that they are satisfied by the basic quantities alone and by the total 
quantities (sum of the basic plus perturbation). Subtracting the equations 
for the basic flow from the total equations and neglecting products of 
perturbation quantities yields 

V' od +U S26! +ad' = FQu', (9) 
ot dX ox yg 


Qa + Ud Gd +vigd+ew=-O, (10) 
a pet Ox oP d+ oP P 





Equation (10) is further modified by substituting for the third 


term 
rag .- —- oo ft (11) 
24 2p ax dp 
yielding, 
ad +KAXb -2¢d K+ 6w'=-Q. (12) 
Bier XP OX oP P 


A diagnostic equation for the vertical velocity is obtained by 


differentiating (9) with respect to p, taking the Laplacian of (12), and 


subtracting, resulting in 
\ 4 
alu 2d +736 = Hw retu+ Via. 3) 
4og ax apex or” P 


Solutions are now assumed to be of the form 


2 Ae = A (pe, t) sinkx +B lp, +) cos ky , (14) 
d 
~ C Cp, €) sin k x +DCp,t) cochkx, (15) 
Qs @, (Sin hx sh) + @. Coos hee § ), (16) 


Next equations (14) through (16) are substituted into equations (9) 


and (13), and the coefficients of sinkxzang cos kx equated giving 


Plas ae as C. B= 6. oer (17) » 
ot Bik” oP 
Bae. ~ heyAg ee oe (8) 
one Gk*> oP” 
MG exe kc =O g Bic Ak - pass: Lt -~k (@, Sing - @ cos §, ) (19) 
te dp GE a 
DD -6kD =-AgAK AU +45k aA tk (Qeor5, +Q 1G). (20) 
aP~ E.- fo sie Ge” a? on 





Equations (17) through (20) are then solved numerically as an initial 
value problem as outlined by Rosenthal [2] : Given initial values of 
A and B, equations (19) and (20) are solved to obtain C and D. These in 
turn are substituted into (17) and (18) to obtain prognostic values for 
A and B. The computational cycle is continued for the desired time in- 


terval. 





Sie Procedures 
Using finite differencing, the diagnostic equations, (19) and (20), 


can be approximated at the odd levels as 


Cia oC ere 9 (Ue 1k ice aks - -B_,) (19a) 


c= | 


“4(4) 4 (@ sin §~ @ cos §,) = 
— : % . 
Baw KD +D.. F thad kK) g (U,,, -& DA; = Ae) 


+4() (Q, cos +Q,sm5,)= (20a) 


¢=).3,5,4.. ie 
At the even levels C and D are obtained by simple averaging, 
C. = Cong C.- D= DD +D:_ (21) 
t ing b~) : " G+) A 
A = 
a er eee 
Boundary conditions are that C and D at the top and bottom levels 
are zero, implying a vanishing vertical velocity there. 


The prognostic equations, (47) and (j6) are approximated at the even 








levels by 
2: 
dAg HKG B= -(¢&)1 (c ia C. )= G , (17a) 
dt k & TAP 
A] 
2Be =~ he A; (6) OF 22) = CG. (18a) 
ot kK‘ 9 SAP 
Cm 2 ae. o 20 
At the odd levels, A and B are obtained by averaging, 
Beene et oes re oP eee oe, (22) 





me | 
Cx bee reer nce ke, 





A vertical increment of 50 mb is used, affording considerable vertical 
resolution. This yields a 21 point vertical grid with 0 mb as level 1 and 


1000 mb as level 21, as shown in figure l. 


Figure 1 - Finite Difference Grid 








Grid-point Grid-point Equation 
index pressure applied 

i 1 0 

ae 50 vorticity 
3 100 omega 
4 150 vorticity 
5 200 omega 

| 6 250 vorticity 
7 300 omega 
8 350 vorticity 
9 400 omega 
10 450 vorticity 
11 500 omega 

| 12 550 vorticity 
13 600 omega 
14 650 vorticity 
15 700 omega 
16 750 vorticity 
17 800 omega 
18 850 vorticity 
19 900 omega 
20 950 vorticity 
21 1000 


The values of C and D are obtained by a simple relaxation scheme as 





follows: 
oe | = 

ee ae (24) 

ni = eo 
. M4 ; 

F; a ge ae ) (25) 

at jams 
ee as a 19, and ae FS =O, then 
Che B.C. +F: (26) 





DT ED the = 
Gs alee” . 


Equations (17a) and (18a) are solved by an iteration technique as 


follows, 

a teat) =A, (erat LG, Bi catia at 

Bo" (ta06) =B. (*)+a$ Cel 44 6) (44000), Cd 
= 2,4,...20, n= 1,2,.... 


Estimates at time t + At are indicated by the superscript n while 


final estimates of quantities at time t are not superscripted. Initial 


estimates for eo t+at) and ee C 4+Ot) are 
Ctteae) = G, (4), (30) 
9 
(£t+0€) = G. (+); (31) 


The iterative process is continued, recalculating C and D for each 


new value of the superscript, until 


| rial tar)-~ A oe eee a, (32) 
(nei) in) 
| iS - (44+O¢) -~B; (#444) | at, (33) 
When the above conditions are fulfilled, the new prognostic values 
of A and B are then used to calculate corresponding values of C and D. 
The process is continued for the desired time interval. 
The amplitude of the pressure wave and phase angle can then be 


calculated using the relationships, , 


2 = Ra (coskx - ¢_? (34) 


3 j 





eo ee | (35) 


tan j, = A/B, (36) 

The model was programmed for selution by the CDC 1604 computer ef the 
U. S. Naval Postgraduate School. The time step, 4 t, was taken as 30 
Minutes, and the tolerance, tT , was 1 CMe Solutions for the amplitude 
and phase angle were obtained for every 6 — up to 48 hours. 

The heating function was introduced in the layers below 500 mb as a 
function of pressure in a manner suggested by Haltiner [2] . For this 
Study, Oy was s@t equal te zere and 


hoo - ‘ 2 
a Sas al we Feo } (37) 
F sa0 


In (37), ¢o xveprasents the amount of heating being censidered in 
joules/gr sec, p is the particular pressure level under consideration 
below 506 mb, and Peooo is the 500-mb pressure, For the purpose of this 
investigation, G2. Was taken as 4 x 1074 joules/gr sec which is approxi- 
mately equivalent to heating of 400 langleys/day. Haltiner and Wang [3] 
propose that tha exchange ef heat between an air mass and underlying surface 
ty!te: in excess of 800 langleys/day. However, it is felt that the 800+ 
figure represents an extreme, and that the value of 400 langleys/day is more 
representative, 

The model is such that the heating function is constant in time and 
Space for a given problem, while the phase angle and amplitude of the pres- 
Sure wave are set initially and then allowed to develop. 

Various permutations of the model were tested with different phase 


angles for the heating functiong and with different wave lengths. 





In all cases, 6 was set equal to .5 x 1074 cgs units, while <9 was 


9713 eect -l 


equal to 1.6 x l sec “. The wave lengths considered were 2000, 


6000, and 10000 kilometers. 
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4. Results. 

The results are expressed as a comparison between a no heat case and 
heating introduced with various phase angles. 

a. The first model studied had linearized vertical shear with 
UL, = 40 m/sec and C= K() <i) , ¢ = 50, 100...1000. The initial 
value of B at all levels was 97.4 metexin and A at all levels was Zero. 
The phase angle of the pressure wave was therefore zero at all levels. 

With the initial input of heat in phase with the pressure wave, the 
development at the end of a 24-hr period in the lower layers was consider- 
able for all wavelengths, (figure 2). However, in the layer 850-700 mb, a 
decrease in amplitude appears for the very long wave, 10000 km. The long 
wave, 6000 km, showed a similar pattern of development with a minimum of 
development in the same layer, 850-700 mb. 

The movement of the pressure wave, as deduced from the phase angle, 
shows that the heating effect caused an increased movement of the long 
waves in the lower layers, (figure 3). Without heating, the very long 
wave retrogressed slowly at all levels throughout the period. For the 
initial 6-hour period with heating, it retrogressed slowly at all levels, 
then quite rapidly in the lower levels. As seen in figure 3, at the end 
of the 24-hr period, the retrogression in the upper levels is less than 
without heating, while in the lower levels it is considerably more. [In 
24 hours, the long wave progressed approximately one-quarter wavelength 
without heating. With heating, it retrogressed initially in the upper 
layers, while progressing in the lower. After 14 hours it progressed at 
all levels, so that after 24 hours, the net effect is progression. The 


short wave progressed over three-quarters of a wavelength without heating 
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at all levels, and slightly less with heating. 

Choosing the 500-mb level as representative of the upper levels, a 
gradual increase in amplitude with time and heating was noted for all 
wavelengths, (figure 4). 

The same model with the heating function initially 180 degrees out 
of phase with the pressure wave, showed that at the end of 24 hours, the 
two waves increased in amplitude below 750 mb, while all wavelengths 
increased below 900 mb, (figure 5). The 10000-km wave was the only one, 
when compared to the no heat case, which showed a larger amplitude at all 
levels. The 6000-km wave exhibited a marked decrease at all upper levels, 
while the short wave decreased slightly. 

The very long wave retrogressed in all but the lowest layers. The 
short wave increased its progressive movement to approximately one wave- 
length, while the long wave decreased its movement, progressing slowly in 
the lower layers, and retrogressing in the upper layers, (figure 6). 

A look at figure 7 for the development at 500-mb clarifies the case 
of the very long wave. It too decreased in amplitude until after approxi- 
mately 18 hours had elapsed. This was true for all other levels. 

For initial phase angles of the heating function of 90 and 270 degrees, 
the same marked increase in amplitude in the lower levels is evident, (fig- 
ures 8 and 9). In the upper levels, for an initial angle of 90 degrees, 
the short wave decreased nearly uniformly, the long waves increased. For 
an initial angle of 270 degrees, the opposite appeared. 

b. A second variation of the basic model included a linear 
vertical shear to 500 mb, and then a parabolic variation of the zonal 
wind to the top of the atmosphere. In this case, Sr equaled 20 m/sec, 


and CZ 


ose WAS 30 m/sec. A further modification was to allow for tilt of 
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the pressure wave by varying the phase angle in accordance with the 
e\ee 
relationship, phase angle equals -7i/a. (Fe) ee, Si 1Q00, 

In all cases, this model showed extreme development below 900 mb for 
all wavelengths, regardless of the initial heating angle. 

When the heating phase angle was zero, the 10000-km wave developed 
very slightly at upper levels, and decreased in the layer 650-800 mb, (fig- 
ure 10). The 6000-km wave developed well above 700 mb, but with a minimum 
of development at 800 mb. The short wave evinced a nearly uniform increase 
in development at all levels above 800 mb. 

For heating with an initial 180-degree phase angle, the 2000-km wave 
showed slight additional development above 700 mb except for a slight de- 
crease in the layer 300-200 mb, and the 10000-km case showed increased 
development throughout, (figure 11). Again, an examination of the 500- 
mb level reveals that the very long wave did decrease in amplitude until 
beyond 12 hours, (figure 12). 

A 270-degree heating phase angle produced increased amplification for 
all levels for the short wave, while both long waves decreased above 800 
mb, (figure 13). 

A 90-degree heating phase angle resulted in an increased amplitude of 
the 6000-km wave at all levels, (figure 14). The 10000-km wave showed a 
decreased amplitude from 900-800 mb, and increased above 800 mb, The 
short wave had decreased development at all levels above 900 mb. 

c. A number of other cases were studied, however, the results 


were comparable to those reported, and therefore are not shown. 
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De Conclusions. 

The restrictions imposed upon the basic model lead to certain 
meteorological improbabilities, i.e., excessive speed of movement in 
certain cases and over amplification after twenty-four hours. However, 
it is felt that for short time periods, the relative development in the 


layers above 900 mb is plausible even on a quantitative basis. 
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APPENDIX I 
Figures 2 through 14 referred to in text are included here for 


easy reference. 


15 











_ os 


es 
ae 


_ 
0% evad balwlonl sve 3x85 At OD tevueteay al devote $ ae ate ; 


iw 
,e2essie? 


“ot 
wih 2» 


CO} 








ho & 
G 


T{ OCO9 = 


ry 


a 


"7 “AVM 7MCOCG 


WHHS VU 





- 
—s 





ARE 


2 aNnora 
SHUNOH 12 LY SYAIDW NI santTtanay 


Oo — 0S2 002 








RUS 9 | 








Alay 












{uy 


HIT dIiOS = ONTLLVH ON 


0 


S 


L 


6) 


Ot OS 


Seamed ay lotic bya 


Ooo t 


006 


008 


ool 


OCE 


O0¢e 


Oot 


0 


aH dH Oe wm 


rH i 


A MMnNND we 


‘Art: 


ON ONS mths patie are —= 


a 





Ree ee Ce Ee ie Oe 


SUHNOH tle LV 





: —— 


il 
q 
= 


SAYIA WH OCC? 
‘IAW WY ICCOO = ie YM HO CoGOt |= "I 
SET a ONTLVSH HLIA 

SGNIT ALIOCS = ONILVEH ON 


rir 


Se te ie Seam Nee i ee ah gts 2 tele aa oa le ay 2: GREE BO eR Sah ac lon Tes RSs Thee tes of ik eed 


© qynold 
SHANI NI AIONV ASVHd 


09 





em a EE tree Coe enti” ena ot 


oa 


O= SI9NV SSVEd ONTLVaH 


O01: 


6) 


an 


Me MRNNOAH R&R 


SHAAHAaAM NH 


A-2 





THART GN 00S 


C—- yr! 
AVA Wy 


Sexil GaHSvd = Oli lv 


Suir] GLiIOS 


V dunold 
SUNOH NI GNIb 


UKs 





0Ot 


O0¢ 


OCE 


CCty 


OCS 


OCS 


ool 


SMHMMmMN 


b+ 


SaaS A 


— Oo ss 





G aynoiad 


SHNOH tie LY SYHIGA NI ednciTdWv 


"MAVM FOIOOO? 
“AAVM WY CCO9 = "I “AVM WYOOCO 
\ AHS = DELLA 
SANIT GIIOS = ONILVEH ON 





- 


CS 00d ost 


OcOt 


008 


ool 


CO9 


00S 


Och 


CO0d 


Oot 


He mer OGMm YN 


UMOEMMNNOM HY 


= 
>= 





Q qunold 
SHNOH qe LY sqayuodd NI WIONV ASVHd 


Oot 


aaAyM WA COCOt 
IT GaHS¥d = ONILVEH HLIM 
SGNIT GIICS = ONILVGH ON 





~- O8t 
COO lt 


006 
008 


ool 


00S 
4 - 
a . 


C02 


7G 
JL= GIONY ESVHd ONTLVEH 


COL: 


fou fom |es) (dp vtep weno emes | 


bH 


Smrirm~A he Oat Ge” 


~A~S 





THATT GW 00S 





\ 
\ 


NST 


\ 
| \ 
LIN \ 


ray 


Td Gensyd =] oven Bae 
S4NTT GLIOS = ONILVGH ON 


Z aunold 
SHNOH NI AWIL 





Ve“ 


LL= BION’ ASVEA ONLLVIR 


00d 


OE 


© 
CG) 
at 


© 
© 
uN 


© 
O 
oO 


ool 


r 


=MWAWMMmuNn 
A-6 


H = 


Gate DAM 





05h 


CON 


Q aunola 
SUNOH NZ LV. SYRIA NI santriay 


OSS OOE cSa 002 


AVM WiCOO0d 


“ZAVA SY CCO9 = 7 rayn *DICOCO 


pereree s 


PVets) 16 ae har RCS) 
SaNIT GITOS = ONILVSH ON 


cot OS 





Ce ALIA Geena) Gti ga 


O00t 


006 


008 


ool 


CO9 


00S 


OOF 


O€ 


© 


00d 


001 





O 


Hz =H aH Ot om UW 


AMM Dw A 








SHNOH 2 


AVA }AICCO?e 

AVIA AYCOCO 
Wek ovl = UnidVeH Hbir 
Sahl] G1i0S = ONLLVaH ON 





eeu L 


6 sunold 
LY. SYSISW NI SALI 


Oo 
=~ PLTONV as Vrid 












Ure es 


Meck 








0 


He =A MyM 


MOM DG w& 


2 





OL sunoLa 


SHNOH Ye LV. SYA NI eantlidny 


0S h 


WliC00¢ 

*TICOQCO L 
wilVin Eibis 

= ONTLVGH ON 


cS2 002 


€ 





COOL 
006 
008 


ool 


C09 


oon 
OCE 


COd 


\ OOt 
a) Ol 
4 a | 


O TSS Cnd Orinbya O 


rH SHAH MOG mY 


WD me A 








WNL U noy (a 


SanlT dLios 


SHNOH 7¢ GV SYSLHW NI sani Tdny 


OO 


I 
nt 
—] 


WOOO? 
7DICOOG Sait 
BiLVen 


i 


= ONITLVGH ON ° 


CSc 





LL ayundid 


00¢ 


oSt 


colt 


ml See 
» 


al GS pegs tio Vid 


Bish 


= 


Von 


006 


008 


ool 


009 


00S 


GOR 


OO& 


CO0¢d 


OOt 


0 


+ a on ee eo 8 oo ep) 


AMM De Ww 


A-l0 


=} See See See 


neem oe en 








THATT GW 00S 






SIAVM WOM OCC’ i 
77 “SAYM PO COOOT! = “7 
Geet GekSVd = ONICVaH HELIA 
SHNTT GIIOS = DNILVEH ON 


cl Sunol 
SUNOH NI GWId 


2 





et 


12> AIONVY dsVid ONTLVaGHR 


OOt 


O00¢ 


Oe 


O 
© 
at 
SwMaHmman 
lo 


rt 


© 

Cc) 

UN 
on 


O 

O 

We) 
GeaanHaDA 


OOL 





CL gyno 
SHNOH 42 LV SUSISN NI SanuiIdyy 


CSE OoOg€ CS2 002d 0S t col cS 





008 


ool 


CO9 


00S 


© 


Gh 


OCe 


00¢ 


O01 





_ HAYM WOOO? 
“AAWA WY 0CO9 = 71 “HAVM HDIOCOOL 
N Gens vd = Ue lwv4H Holi 
SANIT GIIOS = ONILVGH ON 





HA BHAA H MA MW 


a MmMunnDo ew 


het 





ae 





Li 


PL synorIg 


SUNOH To LV SHAM NI SanLlIdvv 


ral 
Aa 
q + Y 
SANIT GIIOS = ONILVAH ON 


- 





CS2 


COOL 


008 


ool 


CO9 


00S 


CCh 


OCe 


C02 


OOt 





tH oe oe ee 2 Oe ep 


we ANN PD & A 


A+13 























